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I. INTRODUCTION 



As has recently been discussed in Ref. the reaction of negative muon capture by the 
hehon, 

/i" + ^He ^ z/^ + (1.1) 

is at present potentially able to provide us with the induced pseudoscalar coupling gp with 
nearly the same precision as from capture by a free proton. Indeed, accurate three-nucleon 
wavefunctions are now available and the uncertainties due to the description of the nuclear 
states can be reduced to a minimum. 

However, the calculations reported in Refs. for other weak reactions on light nuclei 
show that besides the one-nucleon contribution, meson exchange current (MEC) effects 
should be taken into account. That this is so can also be seen from the analysis performed 
in Ref. [0, where Congleton and Fearing compared the results obtained using the elementary 
particle model (EPM) to the standard calculations in the impulse approximation (lA): the 
effective magnetic couplings Gp and Ga are 10% smaller in the lA (see Table |). Table |I| 
shows the contributions of the current components to the effective couplings, and it can be 
seen that j a makes an important contribution to Gp and Ga- 

Here we continue the study of the characteristics of reaction ( p. . 1|) and consider the weak 
MEC effects. For the operator of the weak axial nuclear MEC we adopt the one published 
recently in Ref. 0], which was applied earlier to negative muon capture by deuterons in 
Ref. 0. The weak axial nuclear MEC operator satisfies the nuclear continuity equation 
(PCAC) up to the order considered, which is 1/M^ (M is the nucleon mass). The space-like 
component of the operator contains both static and velocity dependent parts. We take into 
account fully both the static part and the terms linear in external momenta from the velocity 
dependent part. 

Further, we use for the vector part of the weak MEC operator the standard isovector 
currents well known from the electromagnetic processes f^-^. We employ non-relativistic 
operators satisfying CVC at order 1/M. 



A brief account of the work has been recently reported in Ref. ||T0[- Here we present 
the full results. In section y, we discuss the formalism employed, section ^ contains 
the numerical analysis of the studied problem and in section |^ we give our conclusions. 
The most important of them is that including the weak axial and vector vr-MEC in the 
microscopic calculation yields results which agree closely with the EPM in their predictions 
for observables. 

In order to make the paper more transparent, we postpone necessary details of the 
formalism and partial results into a series of appendices. 



II. FORMALISM 

To evaluate the effect of meson exchange currents, the ^H ^He weak nuclear cur- 
rent was parameterized by six Q dependent current amplitudes Po, Pi, jq, jadx^ and jx^"* as 
follows. 
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^H; I f I ^He; ) = xL; {po + PiQ ■ ^]xr. 



(2.1a) 



(2.1b) 



In Eqs. ( |2.1| ), and are the initial and final projections of the internal nuclear 
angular momentum, j° and j are the time and space components of the total weak four 
current and the three-momentum transfer to the helion is Q. The current amplitudes are 
not relativistically invariant and are not independent (if we make the assumption that second 
class currents are absent). To see the non-independence note that the EPM parameterizes 
the current using only four form factors whereas there are six current amplitudes. The 
current amplitudes are, however, useful in the formulation because there exists a simple 
correspondence between the non-relativistic current operators and the current amplitude to 
which the matrix element of that operator contributes. 

The total current has a vector and an axial part. The vector part contributes to po,jQ, j 

(2) 

and the axial part to pi,jcr,jx ■ The relationships between these current amplitudes and 
the traditional effective form factors Gy, Gp and Ga are. 



(1) 

X 



Gy 

Gp 
Ga 



Po + jq 

_ .(2) 
jx jx 



Pi 



.(1) , ■ _ ^ -(2) 
Jx ~r Ju 3"^^ ' 



(2.2a) 
(2.2b) 

(2.2c) 



The statistical capture rate Tq and analyzing powers A^,At and for reaction ( |1 . 1|) 
can be found from the effective form factors by using Eqs. (3) (11) (12) and (13) of Ref. [|1|. 

Our model for the nuclear current employs non-relativistic wavefunctions and current 
operators the latter of which are listed in appendix ^ After suffering a multipole decom- 
position of the plane wave factor, the current operators yielded matrix elements with the 
following five forms. 



-'IXL 



(3H II [47r3^,i(xy) ® S^]jjx{Qy/3)I \\ ^Re)/{2J+1) 



1/2 



(2.3a) 



,ZfxL = ('H II [47ry[i{xy) ® Sj:]jUQy/3)I 



d_ 



^He)/(2J+1) 



1/2 



(2.3b) 



XxL = ('H II [Anyt^im ® Sj:]jUQy/3)I 



d_ 



^He)/(2J+ 1 



,1/2 



:2.3c) 



^ZUl' = (3H II [Anyt.ixy) ® S^]k ® — ixiQy/S)I \\ 'Re) /{2J + 1)'/' (2.3d) 



yKlJ 
y^lXL 



(3H II \[4ny[,ixy) ® S^]k ^ ^] ixiQy/m \\ ^He )/(2J + 1)^/2. 



:2.3e) 
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The letter Z takes on values A ... I according to the spin and isospin operators 5*2 and X 
as indicated in Table III. For example, the Gamow- Teller matrix element is 1/4 Fqqq. The 
matrix elements are reduced in the spin-spatial part but not in the isospin part so that for 
example, 



^He II 3r| || ^He) = +1 



(2.4a) 



(=^H II 8x3^ II 3H) = -1. (2.4b) 
The following selection rules were applied to the matrix elements. 

1. Parity. For Z/^j^ this implies that only matrix elements with, / + A = even, are 
significant. 

2. Hermiticity and isospin symmetry. Because ^He and form a good isospin doublet, 
matrix elements Zf^j^ with an anti-hermitian operator are insignificant. This approxi- 
mation is good at the few xlO~^ level since we have for the wavefunctions used here, 
(which have no isospin 3/2 components), 

(^H II 3/3- II ^He) = 0.9998 (2.5) 



and further the probability of the 1 = 3/2 components is 10~^ — 10^^ [|ri|]. For the 
local matrix elements A, C, D, E, F this implies that only J = 1 is significant and for 
B, G, H, I only J = is significant. A numerical example is -E?io/ (a/S-Euo) = 9 x 10~^. 

3. Selection rule for I. If the spin-isospin part of the operator is even(odd) under the 
interchange or particle labels 2 and 3 then only matrix elements for which / is even(odd) 
are significant e.g. A requires / = odd. This selection rule follows from the symmetry 
of the wavefunctions under interchange of particle labels. 

4. A selection rule peculiar to the F type matrix elements is that F^^-^ is zero by inter- 
change of particle labels 1 and 2. 

Further to the selection rules a rationalization of the matrix elements was possible by 
taking into account the long wavelength of the boson mediating the interaction, the low 
P-state probability (< 0.2%) and the ~ 9% D-state probability in the trion wavefunctions. 
With regard to the first point we note that the isovector spectator point-particle range is 
about 1.7 fm. This implies that matrix elements with high values of A are small and that a 
power expansion in y rapidly converges, e.g. the term in jo{Qy/3) contributes at the 3% 
level. The range of the pair coordinate x is of the order of the pion Compton wavelength, 
1.4 fm, as can be seen from Fig. |1|. This is expected because the exchange current operators 
pick out the spin and isospin dependence in the wavefunction which is mostly due to pion 
exchange. Matrix elements with the factor jf{Qx/2) will be small for large values of /. 

Integrals of the densities are given in Table An example of rationalization is to 
consider the ratios -D202/-O000 = 2 x 10~^ and -D022/-D202 = 1 x 10~^. We note that Dqqq 
receives contributions solely from SS overlap. D202 -^022 dominated by SD overlap 
and are thus diminished due to the lower D-state probability by a factor ~ ^(0.09/0.91) = 
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3 X 10~^ which explains the ratio of D202 -^000 • -^022 is further suppressed due to the 
low momentum transfer in the process by the j2{Qy/^) factorQ ~ 3 x 10~^. This explains 
why Dq22 is so much smaller than D\q2- 

For most local currents the leading multipoles plus those within one unit of angular 
momentum were included. This procedure neglects contributions at the 3% level. For the 
largest current (delta excitation) these 3% corrections were included. 

The matrix elements yZi^^^ and yZ[l}/ were neglected since we expect their contributions 
to be small compared to those of x^^al and ^Z^jf . They will be suppressed by about 15% 
due to angular momentum propensity rules. By parity we need an extra factor a; or y to 
combine with V^^. In the case of x this yields / = A = 1 and these elements are suppressed as 



can be seen in Table |^ In the case of y there will be suppression because of the jiiQy/S) 
factor of about 15%. 

The fact that we omit the non-local terms involving Vy implies that our estimate of the 
non-local terms has an inherent uncertainty of about 10-20%. It turns out that the non-local 
currents due to contribute at the 10% level to the exchange currents and so the above 
neglect affects the results for the MEG at the 1-2% level and the results for the total current 
at less than the 0.5% level. 

Each exchange current yields a contribution to one, (or more), of the current amplitudes 
and the general form is, 

j = J dxp{x) f{x) (2.6) 

where p{x) is a nuclear density and f{x) is a 'potential function' which depends on the meson 
which is exchanged, the overall coupling strength and also the momentum dependence of 
the current. The definitions of the various potential functions entering the calculation are 
given in appendix 0. 

The nuclear densities were calculated using wavefunctions found by the coupled rear- 
rangement channel (CRC) method of Kameyama et al. [|T2|. The Faddeev components are 
expressed as sums of gaussians and one can find analytic expressions for the matrix element 
densities as functions of x, the pair coordinate. The densities were evaluated at either 14 
points or in the case of the large matrix elements (first five of Table |V| and the matrix ele- 
ments with a derivative i.e. xZ/xl)^ 24 points with a higher density of points for x < 1.5 fm. 
The reason for calculating a density is to facilitate the calculation of matrix elements with 
different potential functions which necessity arises when meson coupling parameters and 
strong form factor cutoffs A^, Ap and A^ are varied. 

The calculation of the local matrix elements was checked by comparing results for the 
trion isovector magnetic moment, fi^, with those of Friar et al. [Q. In that work, wavefunc- 
tions for the triton were applied resulting from many NN potentials one of which was the 
AV14 potential: only exchange currents arising from vr exchange were considered. With vr 
exchange only, the expression for //^ in our formalism is. 



^Note that one-body currents have a j\{2Qy/3) factor showing that the effective momentum 
transfer for exchange currents is half as much as for single- nucleon currents. 
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where, 



/x.= nm(-)^jl^) n.m. (2.7a) 
= /i^^ + /^r^ (2.7b) 



l(l + ^P_^n)[^]0,l (2.7c) 



and. 



-^^m[^/i]-78^no[f/i]} (2.7e) 

= (-)— /.V{72^ooo[|4 - t/2] + iDUds] + - 54/x.)] 

(c?5 - 5d3/x.)] - -L.El,,[f{d, - 5ds/x^)] 
-VrAifd,] - ^,El,,[fd,]} (2.7f) 

/if = 2Gif^NAfnNN-^j^JJ^^-^-J^[FQQQ[f3 - /2/2] + ^-^202[/2] 

+ f ^000[/3 - /2/3] - 3^/^2'02[/2]} (2.7g) 

where Kp and the proton and neutron anomalous magnetic moments and 

3.706. 

The calculation of Friar et al. ^ used a different potential function for ^1 arising from 
the propagator lSJ'p{Q2)^^F{Q^)^T^NN{Q2)FT,NN{Q'i) which, although an intuitive choice, is 
inconsistent with the equation of continuity. To facilitate comparison we must replace d2 
/25,o?3 — *■ /26 and c?5 — > /27 in Eq. ( |2.71| ) (see appendix for the definitions of functions 



fi) where f25, f26, f27 correspond to the choice of propagator made in Ref. ||^. For the A- 
excitation graph we need to choose f-j^NA = Q^f^/^f-nNN and Gi = 3y/2/10{fip — fin) = 2.00 
so that our coefficients agree. 

The comparison is shown in Table |V]. Our results for the exchange contributions are 
in overall agreement (case (b)) although the comparison is not exact because we have used 
a triton-helion overlap. Also shown are results with bra and ket wave vectors and 
replaced with — r| (case (a)) which corresponds exactly to the calculation of p. These 
results agree well. The agreement gives us confidence in our calculation of local matrix 
elements. We also give results for 8 and 22 channel wavefunctions from the AV14 potential 
with the Tuscon- Melbourne three-body force (cases (c) and (d)) which show the effect of 
adding extra channels to the Faddeev wavefunction. The effect is very small. Although a 
direct comparison of our results with |^ is not possible because a different coupling scheme 
was employed there (our wavefunctions employ Russell-Saunders LS coupling rather than 
jj) our 8 and their 5 channel wavefunctions have similar content and convergence of the 
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binding energy requires 22 and 34 channels respectively. The extra channels change the 
pair, pion and delta contributions by less than 1% in our case whereas the difference in |^ is 
about 9%. This is a manifestation of the faster convergence obtained when using the CRC 
method: the projection of the potential onto partial waves is nearer to being complete for 
a given number of Faddeev component channels for the CRC method as compared to the 
method used in Ref. 

To check the calculation of the non-local matrix elements we used a Peterson type de- 
^]. The device follows from a simple identity for the matrix element of an operator O 



vice 



multiplied by the momentum operator p when the bra and ket vectors are equal. If, 

M = {ip\Op\ij) (2.8) 



then for O hermitian (O^ = O), 



^M = ^^{i;\[0,p]\^) (2.9) 
and for O anti-hermitian [O'' = —O), 

^M=^{ij\ [0,p] 1^). (2.10) 
Applying the above to the operator O = ia^ ■ xt-^ we obtained, 

3 pi _ rpl _|_ /r) tt rpl , o tt rpOll (0 ^^\ 

2 -'^000 " I'-^OOO "T x'-^202 "T ■ \'^-^^) 

The 'tt' superscript indicates that the triton spin-space wavefunctions were used in both 
bra and ket and that (5 = 0. Eq. ( |2.11| ) relates a local matrix element Fqqq to non-local 



matrix elements both with (x-^ooo'a: -^202) ^"^^ without {xPf^l) derivatives. The relation ( 2.11 



was found to hold to within the precisionQof the calculation and to 0.4% when the ket was 
replaced by a ^He wavefunction. This deviation is first order in small differences between 
the ^He and wavefunctions whereas the deviation given in equation ( |2.5| ) is second order 
which explains why the latter is smaller. 

A further check of the calculation was made by calculating the Gamow- Teller matrix ele- 
ment, M(GT), which measures the axial current at zero momentum transfer and governs the 
/9-decay of the triton. Our results are compared to those of Adam et al. P] in Table The 



calculation reported in Ref. ||^ used axial currents which, except for the A-excitation cur- 
rents, are equal to those used here and applied various wavefunctions and meson parameters 
(the values given in the table are for the Paris potential 10). For the purposes of Table |V| 
we have adopted the meson parameters given in Eq. (3.2) of Ref. 0]. However, because our 
hadronic form factors are monopole FiQ"^) = [(A^ — m^)/(A^ + Q^)]" with n = 1 and those 
used by Adam et al. had n = |, we applied equivalent values for A found by equating the 
slope of F{Q'^) with respect to at the on shell point = —m?. This procedure yielded, 



^The expansion coefficients for the basis functions are know to a finite precision. By noting the 
largest contribution to the matrix element one can calculate the precision of the matrix element. 
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A^(monopole) = 1.69 GeV ^ K{n = \) = l.2 GeV 



(2.12a) 



Ap(monopole) = 2.72 GeV ^ Kp{n = \) = 2.0 GeV. (2.12b) 

Taking into account the different wavefunctions used the results agree quite well. 

To check the calculation of the A-excitation currents we compared our results for M(GT) 



to those of Carlson et al. [|T5[. That calculation used wavefunctions found from the AV14 
NN potential and the Urbana three-body force. Our wavefunction should therefore be 
only slightly different since the three-body force has only a small effect. To make the 
coupling coefficients the same we set /^^^/47r = 0.079, /^atv/^tt = 0.2275, ^'^(O) = —1.262, 
q'pNn/^'^ = 0.5, = 6.6 and Gi = 3.06. The pion and rho cutoffs were set to A^^ = 
0.90 GeV, Ap = 1.35 GeV which makes our potential functions agree exactly since Carlson 
et al. used monopole form factors. With these parameters we found that M(GT) received 
contributions +0.052 and -0.022 from the currents ANP4 and ANP5 respectively which 



agree well with the values reported in Table I of Ref. |15 . 



III. RESULTS 

We precede the presentation of results with a discussion of our choices of parameters. 
The ranges used for the parameters which define the exchange currents are listed in Table 
VTTl Our aim is to make a realistic estimate of the theoretical uncertainty and so we choose 



large ranges of reasonable values. 

We took the ratio of the rho anomalous to normal coupling, ny, to vary between 3.7 
and 6.6. In the context of the hard-pion model, which combines vector meson dominance 
(VMD) with chiral symmetry, the value 3.7 is consistent as this is the VMD value. However, 
the Bonn meson exchange force model OBEPR |]16[ requires Ky = 6.6. A similar criterion 



was applied to find the range of the pNN coupling where the hard-pion model requires the 
KSRF value g^j^j^/in = 0.70 whereas OBEPR favours a stronger coupling of g'^p^^/An = 
0.95. The value of the pNA coupling, Gi, is taken from experimental data for the M1/E2 
multipole ratio for photoproductions of pions |T8|. The parameters i> and 7] reflect the off- 
shell ambiguity in the pion exchange potential. For consistency with the AV14 potential, 
which is a static potential, their values should be taken to be 1/2. Our exchange current 
APSPV expresses the difference between what would be obtained using pseudoscalar and 
pseudovector vrA^A^ coupling. Pseudoscalar coupling corresponds to A = 1 and pseudovector 
coupling to A = 0. It is not possible to achieve exact consistency of this current with the 
NN potential but the value A = 1/2 is the most appropriate. We varied A between and 1. 
The value of f^j^j^/Air should be taken as 0.081 to be consistent with the AV14 potential. 
However, the range 0.075 - 0.081 was used to estimate the uncertainty due to /t^nn where 
the low value comes from a recent Nijmegen analysis of NN scattering data |jl9| . 



The most influential parameters are the strong form factor cutoffs A^,Ap,Aa and the 
ttA^A coupling f-^NA- The strong form factors affect the potential functions fi{x) entering 
the reduced matrix elements and in general reduce their strength at small x i.e. x <~ 1/A. 
The value Am = t^m where M = 7r,p or ai corresponds to complete cutoff. For example 
the current for A-excitation due to p exchange, ANP5, would be zero if the choice Ap = nip 
was made. We used the range A^^ = 1.0 — 1.5 GeV which corresponds to a central value 
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of 1.2 GeV with a 20% variation in 1/A,r either side. The values of Ap and were then 
found by demanding consistency between the exchange current and the NN potential used 
to construct the wavefunction. The exchanges of vr, p, ai yield isospin dependent spin-spin 
and tensor interactions (V5 and Vr) and it is these two components of the AV14 NN force 
which were fitted. The long range parts of Vs and Vr are fitted exactly provided the same 
pion coupling f-^NN and pion mass are used. The short range part is affected by the strong 
form factors and the desired effect of the form factors is that the strength of the potential 
at short distance is reduced. Monopole form factors achieve this effect for Vr but for V5 
an undesirable change in sign at small distance occurs. This can be traced back to the 
second derivative of the potential function arising from vr, p or ai exchange. To illustrate 
this consider the vr OBEP potential (see appendix |^ for further details). 

F' 

V^ ~ F; - ^ (3.1a) 

V-r^F:+2^ (3.1b) 

where F' = dF/dx^. Fig. | shows , V^ , F^ and for A^ ^ cx) and A^ = 1.2 GeV. The 
function F^ changes sign at Xt^ ~ 0.16 and we also have, 

x^fI ^ Flasx^O. (3.2) 

This explains why the combination of F^ and F^ in Vr does not change sign at small 
X despite containing F^ . The combination (Vr — Vs) eliminates F^ . In order not to be 
sensitive to the short range part of F^ we fitted Ap, A^ to Vr and (Vr — Vs) where the latter 
combination eliminates F^ . A posteriori we noticed that the largest exchange currents are 
ANP4 for jo- and VPl for These currents have the form, 

Aj,[ANP4] ~ J dx[F^,,{x) - lDl,,{x)] x jo{Qx/2) V^{x) (3.3a) 

Aj(i)[VPl] ~ I dxDl,,{x) X 3,{Qx/2) x^[V^{x) - VJ(x)] (3.3b) 

The above observation shows that the parts of the potential which are most important to fit 
well are Vf and — in the ranges where (F2Q2 — \D\Q2)jo{Qx/2) and Dqqo ji(Qa;/2)xT 
are appreciable. The above considerations led to the following practical procedure. We 
minimized the function fP'^{Ap, A^) where, 

fP^ = (/^ - /0)2 + (/^^ _ J0^)2 (3.4a) 

and, 

POO 

It = / dxxPV^''''''{x) (3.4b) 
Jo 

4= / dxx'^V^'^^^x) (3.4c) 

JO 
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/■oo 

Its = dxx" {V^''''''{x) - ^^^^(x)) (3.4d) 
Jo 



I^s = I dxx'^ {Vf^\x) - V^^^\x)). (3.4e) 







The choice (p, g) = (3, 3) makes the integrand in It peak at x = 1.6 fm and the integrand 
in Its peak at 1.9 fm and gives the best matching of shape to the nuclear densities for ANP4 
and VPl. Examples of fit values for Ap, are given in Table |V111| . We see that Ap is smaller 
when is larger as one would expect. The value of A^ is poorly constrained and apparently 
may vary between and Ap. 

To summarize, Ap and A^ were chosen to match the isospin dependent spin-spin and 
tensor potentials used to construct the wavefunction. They are functions of A^, fnNN,gpNN 
and Kv and functionals of V^^^^^ and Vg"^^'^. This somewhat artificial procedure would 
be unnecessary if the wavefunctions were constructed using OBEP potentials found from 
the same lagrangian as the exchange currents were derived from. In that case A^,Ap and 
Aa would be fixed by the deuteron properties and NN scattering data. An alternative 
improvement which could be made is to choose a different type of form factor, the idea 
being to match the shape of V^^^'^ and V"^^^^ very closely. 

There are four estimates of the vriVA coupling f^NA- The simplest constituent quark 
model yields /t^na = (^V^/5fnNN and hence f^N^/^'^ = 0-23. Dispersion theory yields 
/^^A/4vr = 0.29 g. The A width of 120 MeV implies f^Nj4:7i = 0.35 The highest 



estimate is /^tva/^tt = 0.36 which is the value implied by relation f^NA = "i/ ^J'^f-wNN coming 
from the Skyrme-Soliton model with l/N^. corrections [^. We note that this model also 
yields /^attv/^^ ~ 0.080 and qa = —1.28 which are in good agreement with experiment. 

We take the range f^N^/4:7i = 0.23 - 0.36 reflecting the various models: there is a large 
uncertainty in the value of this parameter. 

Our calculation does not treat the effect of the A-isobar explicitly i.e. there are no ANN 
components in the wavefunction. We were therefore unable to take into account direct 
coupling to A-isobars present in the nuclei or the indirect effect of the A-isobars on the 
coupling of the nucleons. These processes have been considered for the axial current at 
zero momentum transfer by Adam et al. [^. The direct coupling contributes +0.030 to the 
Gamow- Teller matrix element, M(GT), and the indirect effect -0.021 P^. These two effects 



compensate each other although their sum, +0.009, is not insignificant. 

A further deficiency of the calculation is the static approximation for the propagator of 
a AN pair. Hajduk et al. have shown that this approximation is not valid and leads to an 



overestimate of the A-excitation current by a factor of 1.9 for both M(GT) and /ly, ||2^ , 
Table 1] . The reason for the overestimate is that Ma — M is a poor estimate for the energy 
of the NNA configuration minus the energy of the NNN configuration. This is due to the 
fact that A-excitation occurs only in total orbital angular momentum 2 channels: for these 
channels the kinetic energy is large which contradicts the zero value given to it in the static 
approximation. 

To improve the static approximation we replaced — M by Ma — M + (T) where (T) 
is some average excess kinetic energy. The value (T) = 110 MeV was used, chosen so that 



the contribution of static approximation A-excitation to M(GT) found in Ref. [Q, 0.055 
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is converted to the exact value of 0.031 plus the contribution from the ANN components, 
0.009, i.e. 



X 0.055 = 0.040 (3.5) 



Ma - M + (T) 



This procedure corrects the contribution of the A-isobar at zero momentum transfer but 
the Q dependence of our A contribution is not correct. According to Fig. 3 of p^, we 



overestimate the A-excitation current for Q < 4 fm~^ although the error is small at Q 
0.52 fm~^ which is the momentum transfer for nuclear muon capture by ^He. 

Finally we list the values used for other constants entering the calculation. We have 
used the three-momentum transfer Q = 103.22 MeV, the energy transfer given by lepton 
kinematics Qo = 2.44 MeV, = 138.03 MeV, Mn = 939 MeV, Ma = 1232 MeV, rup = 
770 MeV, = 105.66 MeV, = 92 MeV, ^a(O) = -1.257 ± 0.003, gv = 0.974 ± 0.001, 
qm = 3.576±0.001, gA = — 1.236±0.003 (the values for gy, gu and gA are at = — 0.954m^). 

Besides the uncertainty in the parameters, we need to take into account the fact that there 
are many possible realistic NN potentials. Wavefunctions derived from these potentials will 
yield different matrix elements according to the relative strengths of the different parts of the 
potential e.g. tensor interaction. To take this into account properly requires calculating the 
observables using wavefunctions derived from all the different potentials. We were not able 
to do this but we did estimate the 'model uncertainty' by varying the size of the dominant 
matrix element densities by a constant factor i.e. independent of x. The size of this factor 
was determined as described below. 

The one-body currents are dominated by [1]° and [cr]"'^. The variation in the value of 
[1]° can be neglected since all models will agree at Q = Qand will have the same deviation 
from that value at low Q provided the isovector radius is reasonable. This condition will be 
satisfied for wavefunctions derived from realistic potentials which possess the correct binding 
energy because of scaling . 



The above assertion was tested by calculating [1]° at Q = 103 MeV with the AV14 
and AV14-I-3BF 8 channel wavefunctions which have different binding energies. The bessel 
functions entering the matrix element can be expanded at small Q as, 

jo(2Qy/3) = l-gV/54 + ... (3.6) 

and so the change. A, in [1]° as Q changes from zero to a small finite value is < y'^ > /54: 
which scales as 1/Eb- If we use the mean ^He/^H binding energy and the value [1]° = 0.851 
for the AV14+3BF wavefunctions then scaling implies, 

A{AVU) = 0.149 X + = 0.163. (3.7) 

The above scaling argument predicts [1]'^ = 0.837 for the AV14 wavefunctions which is 
close to the calculated value of 0.839: this result supports the argument that the model 
dependence in the value of [1]'' is small. 



^The value is one minus a correction of the order of few x 10 ^ due to isospin symmetry breaking 
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For [cr]*^'^, however, differences do exist at Q = 0. The impulse approximation contri- 
butions fiy" and M(GT)^^ are both proportional to [cr]^'^ which allowed us to gauge the 
variation in [cr]^'^ due to the use of various wavefunctions. We used values for fiy" reported 
in [| for the Reid soft core (RSC), RSC+TM, RSC+Brazilian(BR), AV14, AV14+TM and 
AV14+BR potentials. Schiavilla et al. report values for /i^ for the AV14+Urbana VII 
and Urbana+Urbana VII potentials. The value of [crf'^ for the Paris potential was taken 
from the calculation of M(GT) reported in Calculations of M(GT) have also been per- 
formed for the Paris, super-soft core, AV14 and RSC potentials as reported in Ref. p^. The 



range for [cr]°'^ from these sources is (-0.913) - (-0.932) and so in our calculation we varied 
[cr]°'^ by a factor ranging from 0.99 to 1.01. We did not include the large values of 0.937 
and 0.943 reported in Refs. for Bonn-type potentials in this analysis because of their 

very different nature. Bonn-type potentials will yield a peculiar balance between one-body 
and two-body currents: basically less two-body and more one-body because of their weak 
tensor force. 

For the exchange currents the dominant matrix elements are -Dqoo' -^ooo' -^202 ^202 
which enter into the A-excitation current and the pion pair (with PS coupling) or contact 
(with PV coupling) term. By studying the 34 channel entries of Tables IV and II of Ref. @ 
we arrived at a common variation factor of 0.917 - 1.083 for -D202 ^"^^ ^202 ^^'^ 0.941 - 1.059 
for D^oo and F^^^. 

These variations due to 'model dependence' made significant contributions to the total 
uncertainty quoted as did the variations in fj^NA and A^r. The uncertainties in the observables 
were found by a Monte- Carlo analysis where the parameters were chosen randomly from their 
ranges with a fiat probability distribution^. It was checked that the probability distribution 
for the observables was close to a normal distribution with the same mean and variance: 
one expects this because of the central limit theorem. 



The final results are shown in Table |X| and the uncertainties listed are one standard 
deviation. Also shown are results for the 8 channel wavefunction: there is little change in 
going from 8 to 22 channels for the wavefunctions used here although we did notice significant 
fractional changes in the value of the smaller matrix elements, in particular those dominated 
by total orbital angular momentum 1 states. We also include results found from the AV14 
wavefunctions. The large difference in the rate is due to scaling rather than being a direct 
consequence of the three-body force. If the one-body currents are taken from the 22 channel 
AV14-I-3BF wavefunction and different wavefunctions used for the two-body currents, then 
the results shown in table ^ are obtained. These results are very similar to each other 
showing that the exchange currents are less sensitive to scaling than the one-body currents 
and that both one-body and two-body currents are insensitive to the increase in the number 
of channels in the CRC wavefunction. 



In table ^ we list the contributions of each current to the current amplitudes and 
effective form factors. The largest exchange current corrections are from the axial delta- 
excitation currents (ANP4-6) and the vector yr-pair (contact) term (VPl). The agreement 
between the microscopic calculation and the EPM is very good except for j^^^ and jj^^ which 



differ by 7% and 11% respectively. Table PCII] shows the separate contributions of local and 



^In the case of experimental uncertainties the range was taken as ib\/3cr 
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non-local currents. 

The dependencies of the observables on the nucleon pseudoscalar coupling gp are shown 
in Fig. ^ and the sensitivies are listed in Table pClll| . The dependence on gp is similar to that 



on the trion pseudoscalar coupling Fp found in Ref. ||T[. The rate is less sensitive to than 
Ay which in turn is less sensitive than either At or Aj\. The curves shown in Fig. |^ are well 
reproduced by parameterizing the effective couplings as follows and using Eqs. (3), (11), (12) 
and (13) of Ref. [§. 

Gv = 0.835 (3.8a) 
Gp = 0.231 + 0.370r (3.8b) 
Ga = 1.300 (3.8c) 

where r is the ratio of gp to its PCAC value. 

Our result for fiy is —2.52 ± 0.03 which agrees with the experimental value of -2.55. 
Our result for M(GT) is 0.977 ± 0.013 which also agrees with the experimental value of 
0.961 ± 0.003. 



IV. CONCLUSIONS 

The calculation presented here used nucleonic and mesonic degrees of freedom to describe 
the charge changing weak nuclear current of the trion system at low momentum transfer. 
The two-nucleon component of the current is given by the vr-MEC obtained from the hard- 
pion lagrangian of the NAnpAi system. The nuclear system is described by wavefunctions 
derived by the coupled rearrangement channel method from the AV14 NN potential with 
Tuscon- Melbourne three-body force. 

We first checked our numerics by calculating the trion isovector magnetic moment and 
the Gamow- Teller matrix element. The results of Table |V| and Table |V| agree well with the 
results of Refs. and 0, respectively. 

Our analysis of the observables for reaction (|1 . 1|) shows that the main contribution 
comes from the space-like component of the current. 

The potential vr-MEC, connected with the nuclear OPEP via the nuclear continuity 
equation, is relatively well defined because the parameters of the OPEP are well known: 
every realistic potential should respect them. This statement is weakened by the fact that 
the AV14 potential is not of the OBEP type and the needed value of the cutoff can be 
extracted only approximately (Table |VI1]) . The axial part of the potential MEG (entries 2-7 
of Table [XID contributes significantly to Gp, while its contribution to Ga is only a minor 
one, because of a destructive interference of the individual contributions. The vector part 



of the potential MEG (lines 14-15 of Table pClD contribute both to Ga and Gp, with the 



prevailing contribution from the pair term VPl. 

There are two sets of non-potential currents. One set is present only in the weak axial 
MEG and arises due to the interaction of the weak axial current with the TipAi system 
(currents ANPl-3). The contribution of these currents to the observables is only a minor 
one (lines 8-10 of Table [XH) . The other set of non-potential currents is formed by the A- 



excitation currents ANP4-6 and VNPl-2. The axial currents ANP4-6 contribute to Ga 



considerably (lines 11-13 of Table PQ) . This set of currents is much more model dependent 
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and it is mainly responsible for the uncertainty of the calculation. It follows from the analysis 
of section III that the contribution from these currents is slightly overestimated due to the 
static limit in which these current are considered here. 

One of our main goals is the analysis of the possibility of extraction oi gp. We have 
found that the rate is rather insensitive to it but the spin observables offer the opportunity 
of measuring gp precisely (see Fig. |]). However, these kind of experiments are also more 
difficult to perform. 



We presented the final results for the observables in Table It can be seen that the 
microscopic calculations and the EPM predictions agree very well. In particular, our result 
for the transition rate for reaction (1.1) is. 

To = 1502 ±32 (4.1) 

Our estimate of the uncertainty in the calculation yields an error of ~ 2 % in the value of 
the transition rate. The largest part of the uncertainty in Fq comes from poor knowledge of 
/ttAta • The value of Fq Eq. ( |4.1| ) is in a good agreement with the preliminary results of the 
new precise measurement 



Fo^'P = 1494 ± 4 (4.2) 
Using this value we can compare ( ^.1| ) and ( [4.2| ) and conclude that 

1. the structure of the space-like component of the weak axial vr-MEC is well described 
at low momentum transfer within the framework of the phenomenological hard-pion 
method, 

2. the value of the induced pseudoscalar constant gp is, 

-1^ = 1.05 ±0.19 (4.3) 

yp 

and so the PCAC value of gp is in rough agreement with the data as is the heavy baryon 
chiral perturbation theory value, gp^^/g^^^^ = 1-01 ± 0.02 p9[. By g^^^^ we mean the 
value found from, 



^PCAC(^2) ^ ^^,^(,2) (4.4) 



which yields 5'p^^*"(— 0.954m^) = 8.12. Let us note that the best measurement of ordinary 
muon capture by the proton yields and uncertainty in gp of 42% and combining various 
measurements reduces this to 22% pl[] . 

The 0.19 uncertainty in gp results almost entirely from the 2% theoretical uncertainty in 
the calculation of the rate. If the spin analyzing power was measured and used to infer 
the value of gp/ g^^^'^ then the lower limit on the uncertainty set by theory is 0.02. 
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APPENDIX A: MOMENTUM SPACE CURRENTS 

In this appendix we list the currents used in the calculation. They are written in their 
momentum space representation. Our conventions are that qa < and that the total current 
is the sum of the vector current jy and the axial current Ja, not the vector current minus the 
axial current. The overall sign of the currents are consistent with the dominant one-body 
axial current. 



(Al) 



The one-body currents are listed below and are consistent with the scheme used in [Q. 
Here, the initial and final momenta of nucleon i are written p and p' respectively and the 
Pauli spin matrices for nucleon i are written a. 



Ja 



+9p 



a ■ p 



{P' + P? 



+ 



1 



2M 



4M2 

12 



P 



(ip X p' + p'{p- (?) + p{p' ■ a)^ 



a ■ p 

2M 



P 

8M2 8M2' 



(A2) 



Pa 



T°- r a ■ {p' + p) 



2M 
P 



P 



12 



3p2 

8M2 8M2 



o ■ p 
2M 



(1 



p' 



3p 



12 



8M2 8M2 



(A3) 



fv 



p + p' ia X [p' — p) 



2M 



2M 



+ 9m 



ia X Q Qo / Q iax{p' + p) 



2M 



2M^2M 



+ 



2M 



)]} (A4) 



r 



Pv 



{9v 



{p' — pY ^ ia ■ {p' X p) 



8M2 



4M2 



+ 9m 



{p' — p)2 ^ ia ■ {p' X p) 



4M2 



2M2 



} (A5) 



Now follow the two-body currents which have been labeled APi, ANPi, VPi or VNPi for 
the sake of reference. A(V) stands for axial- vector (vector) and P(NP) stands for potential 
(non-potential) current. The currents are written in terms of non-local momenta P2, -P3 and 
local momenta Q2, Qs defined by. 
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P2=P2+P2 (A6) 
Q2=P2-P2 (A7) 

P3=p's+P3 (A8) 

4=P3-P3 (A9) 

where (pi) is the momentum of nucleoli i in the final (initial) state. The currents appearing 
below are for the pair of particles labeled (23) and the isospin component a e {x, y, z} i.e. we 
have written here j"(23). The total current for muon capture is then j^~*^(12) + j^'~*^(21) + 
f -iy (^23) + f-^y (32) + f-^y (31) + f-'y (13). Given the current j(23) for particles 2 and 3, 
the matrix element of the other five other currents j(32), j(12) . . . are equal to that of j(23) 
which follows from the symmetry of the wavefunctions under interchange of particle labels. 

First we list the weak axial potential currents; 
the sum of the Oi-pole pair term and the PCAC constraint term, 

J^'^-^[AP1+AP7] = (^)' (Q3)i^.V(Q3)(a3 • Q3){[Q + 1^2 X P2K 

+ [P2 + ia2 X Q]i{t2 X rs)"} (AlO) 

where the total current is related to the bare current for this current and currents AP3,AP4 
and APSPV by,g 

Jl = Fa{Q)F^nn{Q)1T'' + F.^iv(g)(l - FAmj^^-TyiiQ ■ jfn. (All) 
the TT-pole pair term plus 7r-pole contact term, 

x{[g' + Q3 + ^^2- (P2 X Q2M 

+ [(Q ■ P2) + 4(4 ■ P2) - 3(4 ■ Ps) + 3ia2 ■ [Q x Q,)]t{T2 x r,r} (A12) 
the TT-retardation term, 

Jr-^[AP3] = - {^[AUQ3)rF^^^{Q,){as ■ 4) 

x{(l - z/)(Q ■ Q3)[Q,r^ + 1^2 X Q3«(r2 x r^Y] 

+ [(1 - V)P2 • 4 + (1 + ^)P3 ■ 4][4^(r2 X + 1<J2 X QgTgl} (Al3) 

the TT-form factor term. 



^In the calculation, the axial form factor Fa{Q) was set to one so that the second term m Eq. 



( All ) yielded nothing. The correct value of Fa{Q) is one minus 0.019. 
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9 



9a 



d 



2MJ 4M^^^^='^rf%^'^^^^^ 



x(a3 ■ Q3)|(l - v)iQ ■ Q3)[Q3rS + X Qst{T2 X 

+ [(1 - r])P, • 4 + (1 + V)P3 ■ QslfeK X r^r + X Q,t^]}, (A14) 

the current which measures the difference between exchange currents derived using pseu- 
doscalar and pseudovector ttN couphng, 



Ja,barej^pgpyj = -(1 - A) 



f 9 



9A 



X ■ 



+ 



\2Mj AM 
{ [[P2 + ta2 X Q]{as ■ Q-,) - Q;{as ■ Ps] 

[Q + id2 X P2](a3 ■ 4) - ia2 X 4(^3 ■ A)] 7-3"} 



H^2 X TsJ 



the ai-pole contact terms, normal plus anomalous, 



Ja,barej^p5j 



9 



2 1 



X 



2My 2gA 
{^3 ■ Qs' 



(Q3) -pTriVAr (Qs) -PpA^Af (Q2 ) 



2M 



{P2 + (1 + /s:y)ia2 X 4)^(^2 X rs)'^ 



(A15) 



(A16) 



x(a3 ■ 



4){i 



8M2 



■[-Ql + ta2 ■ (4 X P2)]}^(t-2 X T^)\ (A17) 



where the relationship between the total and bare current in this case is| 



"?a,bare ^ a, bare \ 

[Q ■ J A - QoPa J 



Here follow the weak axial non-potential currents. They are; 
the OipvT term. 



Ja,barej^^p^j 



f 9 



\2Mj SMqa 



A''f{Q2)AI{Q3)FpNn{Q2)F.Nn{Q3){^3 " 4) 



(A18) 



{ [P2 + (1 + ny)ia2 X 4] [4 ■ (4 - 4) + ^\S^{P2 ■ 4 - ^3 ■ 4j 



M2 



+4 [4 ■ (A - i^2) - (1 + «:y)^t?2 ■ (4 X 4) 

(1 + 2/ty) 



-{P2 ■ 4)[-4 +^'?2 ■ (4 X P2)]]}^(r2 X r3)'^ 



(A19) 



gain, the Fa{Q) = 1 approximation was made so that the second term in Eq. ( |A18 ) yielded 
nothing. 
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f 9 



\2MJ AgA 



AUQ2)AUQs)FpNn{Q2)F^Nn{Q3){^3 ■ Qs] 



(Q2 ■ P2 



(1 + 2/ty) 

8M2 



[-Ql + ^^2■{Q2y<P2)] 



{Q2 ■ P2 



■[{Q2 ■ P2) + (1 + /ty)^CT2 ■ {Q2 X Q3)] K(^2 X T^Y 



(A20) 



where the relationship between total and bare current is given by Eq. ( [A18| ), 
the pvr term, 



j^[ANP2] = - 



2m) J]^^-^^(Q)^F{Q2)AUQ3)FpNN{Q2)F^NN{Q3Kr2 X T^T 



X 



(^3 ■ 4) {^2 + (1 + /ty)^^2 X Q2} 



(A21) 



9 V 



p^[ANP2] = - ( ^i^.7Viv(Q)A^(Q2)A^(g3)i^piViv(Q2)i^.jviv(Q3)^(r2 x t^T 



X(^3 ■ Q3){l + l-Ql + ^^2 ■ (4 X P2)]}, (A22) 



8M2 



the pmr term, 



J^[ANP3] = -g ^^A^(Q)A^(Q2)A^(g3)i^piV7v(Q2)i^.iViv(Q 



:(f?3 ■ 4) {4 ■ {P3 - P2) - (1 + ^^vY^2 ■ {Q2 X 4) 
+ 2k^ 
8M2 



+ ^3 ■ ^3) ^^^,y [-g^ + ■ (Q2 X P2)]Ut2 X ra)^ (A23) 



delta-excitation, tt propagator. 



X( 



:('?3 ■ Q3){4r3" + W2 X 4)^(^2 X T,r}, (A24) 
delta-excitation, p propagator. 



X 



47va\ (1 + Ky) AGig^U 
) 2M2 9(Ma-M)' 

{4 X (^3 X 4)^3 + ¥^2 X (4 X (^3 X 4))^(t2 X Ts)"}, (A25) 



delta-excitation, ai propagator. 



2 



Jl[ANP6] = 9\9' .., _,.. A°(Q3)F.W(Q3){cr3r3° + iza2 x a3z(r2 x r3)»} 



/ttA^A \ 3 2 ^ 

^^^''9(Ma-M)' 

(A26) 
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where, 

as^as + ^^. (A27) 

Finally we list the weak vector currents. They are; 
the TT-pair term, (with PS coupling) or vr-contact term (with PV coupling), 

J^[VP1] = - (i^) Al{Q,)Fl^^{Q,)a,{a, ■ Q^r, x ^T, (A28) 
the pion current, 

J^[VP2] = + (^] 4(^2 • Q2)(^3 • Q3)i{T2 X ra)" 

Xg|3g| [^UQ3)F^NNiQ3) - Al{Q,)F^^M2)] , (A29) 
delta-excitation, tt propagator, 

X(<J3 • 4)«Q X {QsT^ + \i{B2 X g3)^(T2 X Ts)'^}, (A30) 

delta-excitation, p propagator, 

:ig X {Qs X (as x g3)r3" + \ia2 x (Qs x (as x Q^))i{T2 x rg)"}. (A31) 



APPENDIX B: TRANSFORMATION INTO CONFIGURATION 

REPRESENTATION 

Purely local currents of the form /((52, Qz)g{Qz) were transformed into configuration 
representation according to, 

J23(^, y) = exp[ig • (f - f )] j exp[-ig • x] x f{Q - q, q) g{q) (Bla) 

= exp[zQ • (f - f )]/(g - zV„ zV,)^(x) (Bib) 

where, 

f = f2 - fa (B2a) 

y = f - ^ (f2 fs) (B2b) 
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j^--^exp[-iq-x]g{q) (B2c) 

and V^. is a derivative acting only on the function immediately to its right, in this case 

This process yields an operator in configuration space with a plane wave factor, a spin- 
space part made from ^2,^3 and x, an isospin part and a potential function which is a 
function of x =\ x \. As an example consider the vector current vr pair term. 



jviQ2, Q3) = - ( ^] ^{r2 X ra) V2(a3 ■ Q,)AUQs)F^^M3) (B3) 



Jv{x,y) = +f^^^exp[iQ ■ (f - f)]z(r2 x r3)"za2('?3 ■ x)fi{x) (B4) 

The potential function in this case is fi{x), (see Fig. |^), the isospin operator is i{T2 x r^^)"' 
and the plane wave factor is exp[iQ ■ {x/2 — y/S)]. 

Non-local currents {aP2 + l3Pz) f {Q21 Qz)9{Qz) were transformed to, 

J23(^,y) = exp[iQ ■ (f - f)\{f{Q -iVx,iyx)g{x) 

X [a{-2iV:, + iVy) + /?(22Vx. + iVy)] + [a{Q - iV^) + (3iV^]f{Q - iV^, iV^)g{x)} 

(B5) 

An example of this type of term is part of the AP1+AP7 current. 

Ja(4, 4, P2, P3) = - (^] ^Ar^ X T,r X ^(^3 ■ Q,)/Vp{Q,)Fl^^{Q,) 

\ itLt^ J 2M 

(B6) 



Ja{^, y) = -f^NNdA^eM^Q ■ (f - f)]^(^2 X Ts^ X {2(^3 ■ [-2iV^ + iVy] 

+—Qi{(f3 ■ x)fi{x) - x{a3 ■ x)f2{x) + <T3/3(x)} (B7) 

The operator V was written as the sum of a part with a derivative xd/dx and an angular 
part V x/x where, 

V. = 9 A + (B8) 

odx sm Ux ocpx 



APPENDIX C: LEADING CONTRIBUTIONS FROM EACH CURRENT 

Here we write the contributions which were included from each current. The argument of 
the two-body matrix elements are potential functions /«, (defined in appendix 0), multiplied 
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by either jo{Qx/2) which is suppressed for succinctness or ^ 0) which is then written, 
e-g- x5(5oo[/i] means, 

(3H II 12/73 [47r3^°o(a;y) ® So{a2,as)]i'i{r2 x r^y h{x)UQx/2)UQy/3)^ \\ ^He ) (CI) 



and Fl^2[hf2] means, 

(3H II 12/V3[47r3^i\(:ry) ® a3]iT-3-/2(x)j2(gx/2)ji(Qy/3) || ^He). 

The one-body matrix elements [1]°, [cr]^ and [cr]''^ are defined by Eqs. (47-50) of [|T| 
One-body current 



(C2) 



Apo=[gv{l-^)-9M 



4M2 



(C3) 



Q_ 

2M 



-^+9p — (1 
3 m„ ^ 



24M2 



cr 



(C4) 



2M 



3 ^*^2M 



[1]° 



(C5) 



(1) 

a; 



Q_ 

2M 



6M 



(C6) 



lO,l 



1 g 



24M2' 



3 2M 



Q ^ 
.^^3M+^^m 



24M2 



(1 



(C7) 



'1 + 



r^i2,i_ ^ 



9A7nrp + 9p — 1 - 



2ML^ 3M 



24M2 



(C8) 



API + AP7 

Aj. = -gAf^NNF.NNiQ)"^ X {iBUfi] - f^Bl,,[f,] + f^HUfi] + iHUfi] 



AP2+AP6 



Ajf ) = 3Aj. = -gAf'NNF.NNiQ) . 



Q Q 



{Q^-Ql + ml)2Mm^ 

^{~72-^OOo[/2 - 3/3] - |-D202[/2] + |^-Bdoo[/l] " "^^-S202[/l] + "^^^-f^OOO [A] 

+L-^202[/i] + 73^-^202 [/i] -X -Bidi[/i] - Vf^Cioii/i] - VL-f^idi [/l] 
-v/f/m[/i]} (CIO) 
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AP3 



= -^QAf^NN^^NNiQ)^ X {|^fiooo[5/20 " /19] " ^^^202(2/20 " /19] 



+ ^7-^000 [5/20 - /19] + L-f^202[2/20 - /19] + :33„-^202[2/20 " /iq] " t^B^Ollf' 



+ 73^-^101 [/20] - ^^-6101^/20] + ^^-f^?0l[/20] - V'i^-f^l0l[/20] 
~V^^a:-^10l[/20] - VljlOl[f20] + VljlOliho]] 

AP4 the same as for APS with 

ly -277 

/l9 ~^ /2I 
/20 /22 

APSPV 



/20 



(Cll) 



(C12) 



APS 



Ai. = -(1 - \)gAf:NNF.NN{Q)^ X {i-^BUfi] + f.BlM - ^^HUfi] 

~L'^202[/l] - 73j.-^202[/l] + ij;-Sldl[/l] - 273^.-^101 [/l] 
+ V^Ma;-^10l[/l] - V^^H^Olifl] + V^i^-f^lOl[/l] 
+ V^^a;-^10l[/l] + V^i^-^10l[/l] - ViiJlOlifl]} 



(C13) 



^ttNNW)^ X 

^2 



+ r 



v/2L>o'oo[/5/3-/6] + |i^2'o2[/5]; 



where, 



+ ~L-^OOo[/4] + ^^^2021/4] +x BlollfA]] } 

= G3{f ,i^doo[/4] + i^2'02[/4] - VlMolif^l} 



(C14) 
(CIS) 



(C16) 



ANPl 



gA 2Mmp " ^ 2 



X {^^^[^/2i?o'oo[/i5/3 - /i6 - i(Sj'(/i2/3 - /la)] 

+ i-D202[/l5 - |(£^)^/l2] + l^-BoOot^Tr/ie + 5/l8 - /l7 " i(£^)^/ll] 



3 3.-^202 [^^^r/ie + 2/l8 - /l7 " |(£^)^/lij 

^?o\M^./i6 + |/i8 - KSjVn] + ^,5i'o\M/i8] - "f^Bl^lihs]} (C17) 



^^.(2) ^ _f^^^ ^ X {V2/^000[/l2/3 - /is] + |i^^02[/l2]} (C18) 
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+'^3 ^^{+x- ^202 [/is]- ^^Dllllfis] + ^^^loil/is] - ^^^ioi[/; 
Zkv mp ^ ^ 



ANP2 



'-^G3^{^,i?^00[3/l3 - /12 - ^(3/6 - /5)] - ii^2^02[/l2 " ^/s]} (C19) 



^^.^ = -^2^^i^F.iv^(g) X {V2i^^oo[/i2/3 - /i3] + IDIM} 

-^-f^^F^MQ) X {-Li?o'oo[/ii] + ^.i?202[/ii] i??o}[/ii]} (C20) 



771 7r 



(C21) 



ANP3 



Unn Y fm^Y (1 + 2Avy)mpQ^ ^ ^ . 
"i^J ImJ 4MK-gg + Q^)M3.^ooo[/i7-5/i8] 

+^^-^202 [2 /i8 - /17] + |^5?di[/i8] - 73^5wi[/i8] + 



(C22) 



ANP4 



4m, 



9(Ma - M) 



X 



{FoO0[/3-/2/3] + ^F2^O2[/2] 



+ ^^000[/3-/2/3]-^i52'02[/2]} 

Higher order terms for ANP4 



(C23) 



4m, 



_ ^) X {Fo1oo[J2/2] - ;^i=^202 [72(2/2/3 - /a)] - f,FUjo{f2/3 - m 

+ ¥^110[jl/2] + i^ll0[jl(17/2/30 - /a)] + ^F2\o[jo/2] + 275^22l[jo/2] 

+V^^222bo/2] - 472^000b2/2] - !^2'o2 b2 (/2/6 " fs)] - §^5^22 U (/2/3 " fs)] 



(C24) 



24 



-V^^^mb'iM} (C25) 

ANP5 

+^ ^ooo[/7 - /8/3] - ^Aifs] } (C26) 

ANP6 

+ f ^000[^8 + /24 - /23/3] - i^i^2'02[/23]} (C27) 

VPl 

+ 76Cuibo/i] - ViiElM,] - VIeUjoM} (C28) 



VP2 



A£^ = /.V^{;^^ooo[e3/3 - 62] + iDUe^]} (C29) 



VNPl 



AjI^^ - -UnaUnN q^^^^ ^ {-^000 [/3 - /2/3] + ^-^202 [/2] 

+ ^^000[/3 - /2/3] - ^i^^02[/2]} (C30) 

VNP2 

" 9(Ma - M) l^j ]i^xF000[/7 - /8/3J + -F202[/8j 

+f ^ooo[/7 - /8/3] - ^,DUfs]} (C31) 

APPENDIX D: DEFINITION OF POTENTIAL FUNCTIONS 

Here we define the potential functions /j. They are derivatives of root potential functions 
Fi, e.g. /i = —dFi/dxT^, f2 — (fFi/dx^ — 1/x.^dFi/dxT^. We first define the root potential 
functions: 
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1 1 

47rFi = Yo{x^) - R}Yo{xk^) - \rJ{R^ - 1)F(xaJ (Dla) 

47rF2 = \Y{x,) + \R,~^Y{x^^) + -^[i?|yo(a;Aj - Y^{x^)]) (Dlb) 

^^^3 = -^^[R^Y^ixK^) - Fo(xO] + ii?;^y(^^J + ii?^ 2 ^i^^ _ l)4^r,(xAj (Die) 

1 1 

47rF6 = (!!!£)'{r„(i,) - ii|r„(iAj - - l)r(xAj} (Die) 



1 

m^6[Ap, A^]6[Ap,m^] ° m^h[K,mp]h[mp,m^] ° J 



where i?M = {KM/mMY and = exp(— x), Yq{x) = Y{x)/x, Yi{x) = Yq{x){1 + l/x), 
Y2{x) = Yo{x){l + 3/x + 3/x^), Ys{x) = Yo{x){l + 6/x + 15/x^ + 15/x^), b[a, p] ^ - p^, 

X-j^ = TTlj^X, X\^ = A.^X, Xp = UlpX, X\^ = ApX, = TTla^X and 3^Aa-|^ ~ A^j^X. 

We now define the fi in terms of the root potential functions with rules like, 

h = F^[Y^ -^Y2,Y^ yY,; l/mj]. (D2a) 

The rule means that /g is equal to Fq with Yq replaced by Y2 and Y{^) replaced by 
where ^ can be any argument. Further, each term exp(— ax) becomes (a/m^)^ exp(— ax) 
which is indicated by the after the semi-colon. 

fi = F^[Yo ^Yi,Y^ Y; 1/m,] (D2b) 
/2 = F^[Yo ^Y2,Y^ yY,- 1/m^] (D2c) 
h = /i/^. (D2d) 



fA^F^[Y^^YrMm-,] (D2e) 

/5 = F4[yo^l^2;l/m^] (D2f) 

h = A/^:. (D2g) 

/7 = FelFo ^ lo + l^i/y, r ^ y - Fo; l/mj (D2h) 

/s = FefFo ^Y^.Y^ yY,; l/mj (D2i) 

/ii = F7[>o--n;l/m.] (D2j) 

/12 = i^7[i"o ^ Y2; 1/ml] (D2k) 
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/i3 = hilx. (D21) 

/i4 = F,{Yo ^ + Y^/y- 1/m^] (D2m) 

/i5 = F,[Y^ ^ Y^, l/ml] (D2n) 

/i6 - F,[Yo ^ yi/y; 1/mt] (D2o) 

h7^F-j[Y,^Y^-l/ml] (D2p) 

/i8 = (D2q) 

/i9 = F2[Yo ^Y:,,Y^ yY^, 1 1 ml] (D2r) 

/20 = F^[Y^ ^ Fa/l/, ^ 1/m^] (D2s) 

/21 = F^[Y, ^ F3, >^ ^ 1/^-2, yVi ^ yVi; l/m^] (D2t) 

/22 = Fain ^ 1^2/?/, 1^ ^ Z/'n ^ 1^; IM] (D2u) 



/23 = F^[Y^ ^Y^,Y ^ yYv, IjmW (D2v) 

/24 = i^s^o ^ n/y, 1^ ^ l/ml^] (D2w) 

/25 = F2[Yo ^ Fi, r ^ F; 1/m^] (D2x) 

/26 = i^2[l^o ^ 1^2, 1^ ^ l/mll (D2y) 



The functions are given by, 



87re2 = — [eii(m^) - eii(A^) + (i?^ - l)ei2(A^)] (D3) 



Snes = e3i(m^) - 631 (A^) + (R^ - 1)632 (A^) (D4) 



where. 



eii(m) = dt exp{-cx)\jo + ^ji] (D5) 
ei2(m) = dt exp(-cx)[f Jo + f (1 + cx)n] (D6) 

631 (m) = — fdt exp(-cx)[(i + c - 2!^)jo + Qt{^ + 2)ji] (D7) 

e32(m) = dt ^^^Pt^[(ca; - ^^(1 + cx))j, + f (2ca; + 3 + (D8) 

the argument of the bessel functions is Qtx and the m dependence comes from c which is 
given by 
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c{m) = [QH{1 - t) + rri^]^ . (D9) 
The functions di are given by, 



APPENDIX E: CONSTRUCTION OF OBEP 



The exchange of tt, p and ai mesons leads to an isospin dependent NN potential with 
tensor and spin-spin components. The non-relativistic momentum space potentials between 
two nucleons labeled 1 and 2 were taken to be, 



{(Jl-Q2){(J2-Q2)^UQ2) 



mi 



(Ela) 



X 



mi 



—PInn{Q2) { — 



(Elb) 



where. 



3Ca, Ti • T2 



(7i • (72 + 



(CTI ■ (52) (0-2 ■ 



ai 



m. 



(Elc) 



(E2a) 



^7!" 2 /I I \2 / 



2 / 



3 .p.vivv- . -w y2M. 

The configuration representation of these potentials are, (with M — tt, p, Oi), 

V^{x) = CM[V^{x)S^2 + Vi'{x)a, ■ a^] 

where, 

5'i2 = 3(ai • x){a2 ■ x) - ai ■ a2 

and, 

V^{x) = f:- F'Jx, 



V^{x) = -{F';-Fl/xp, 



V^\x)^Fl-F'jxa, 



(E2b) 

(E3) 

(E4) 

(E5a) 
(E5b) 
(E5c) 
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Vs'{x) = -2{F'; + 2Fl/x,j 



(E5d) 
(E5e) 



Vs^^{x) = F:^+2F'Jx,,-3F, 



(E5f) 



with, 



°° gsin(gx) 



(E6) 



The notation Fj^ means d/ cLxmFm- With monopole form factors, FMNN(q) 
m\j)/{A\j + g^) we found that, 

AnV^ix) = (-l)^W^) f^) X [Y,{mMx) - RIM^mx) 
1 

-RI{Rm-1)I^mxY^{Kmx)] 



(A2 



M 



(E7) 



47r V^/(x) = Yo{m^x) - R^Yo{Kx) - Rn^R^ - 1) Y{A^x) 



(E8) 



ArrV^ix) 



-2) ( ^ ) [Fo(mpx) - R^Yo{A,x) - i?p2 _ y(A^a;)] 



(E9) 



4nVs'{x) 



[Yo{ma^x) - RlYo{Aa,x) 



Rai^ 



{Ra,-l){Ra,-S)Y{Aa,x)]. 



(ElO) 



Our exchange potentials for vr and p are the same as those used by others p2| , p!6| . Our 



exchange potential for oi agrees with that of Ref. |32| but not with that of Ref. |3^. The 
difference is that the (cxi ■ Q2)(o"2 ■ Q2) term was left out in |Q. Here we include this term 
because we found that it contributes non-negligibly: V^{x) is entirely due to the term and 
Vg{x) would be ~ 50% smaller at a; = l/A^^ and 50% larger for x ^ l/Aaj without the 
term. 
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FIGURES 



FIG. 1. Examples of nuclear densities. Points were also calculated at x=5.5,6.5 and 7.5 fm but 
are not shown. 

FIG. 2. Isospin dependent spin-spin and tensor potentials from pion exchange and first and 
second derivatives and F^. The solid line is for = 1.2 GeV and the dashed line for A^r = oo. 
The constant Cj^ has been set to one for these graphs. 

FIG. 3. Variation of observables with the nucleon pseudoscalar coupling gp. 

FIG. 4. Potential function for the it pair term of the vector current. The pion cutoff is 
A^ = 1200 Mev 
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TABLES 



TABLE I. Effective couplings in the EPM and lA and their difference. The uncertainty in Gp 
reflects experimental uncertainty only. 

Model Gv Gp Ga 

EPM 0.85 lb 0.01 0.603 ± 0.001 1.29 it 0.01 

lA 0.84 -1% 0.523 -13% 1.19 -8% 



TABLE II. Origin of contributions to effective couplings 





Gv 




Ga 


Dominant 




=5 =J 


JA^ 


Other 


JV 


PA 


PA,jv 



TABLE III. Labeling of reduced matrix elements. The superscript in the isospin part means 



z 




I 


A 


u 


12i{T2 X Ta)- 


B 


^3 


12i(r2 X ra)- 


C 




12i(r2 X ra)" 


D 


[(72 0-3] 1 


12i(r2 X ra)- 


E 


[(72 (8) C73]2 


12i(T2 X Ta)~ 


F 


^3 


12x3- 


a 


[(72 ® (73]o 


I2T3- 


H 


[(72 1^ 0-3] 1 


12x3- 


I 


[(72 0-3] 2 


I2T3- 
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TABLE IV. Sizes of matrix elements. The Bessel function in y has been expanded and only the 
first term kept except in the cases of I?ooo ' -^202 ' -^000 ^202 • The dominant overlap describes 
the total orbital angular momentum of the wavefunction components which contribute most to the 
matrix element. 



rix element 




Integral 




Dominant ovei 




AV14+3BF[8] 


AV14[8] 


AV14+3BF[22] 




-^000 


9.62 


9.57 


9.63 


SS 


^000 


5.78 


5.96 


5.77 


SS 


-^000 


-3.55 


-3.53 


-3.55 


SS 


-^202 


2.03 


1.99 


2.05 


SD 


-^^202 


-0.900 


-0.867 


-0.907 


SD 


-'202 


-0.066 


-0.056 


-0.082 


SD,DD 


-^110 


0.20 


0.20 


0.20 


SD 


-'110 


-0.15 


-0.15 


-0.15 


SD 


-^110 


0.14 


0.14 


0.14 


SD 


^111 


0.057 


0.061 


0.054 


DD 




0.045 


0.053 


0.044 


SS,DD 


^111 


0.035 


0.037 


0.025 


DD 




-0.024 


-0.023 


-0.031 


DD 


-^110 


-0.017 


-0.018 


-0.017 


SS,DD 


-^111 


-0.012 


-0.014 


-0.012 


DD,SP 


-f^iii 


0.004 


0.003 


0.000 


SP 


-'^022 


-0.017 


-0.017 


-0.017 


SD 


-^022 


0.005 


0.005 


0.005 


SD 


-'022 


-0.001 


-0.001 


-0.001 


DD 


xJ^QOQ 


-7.00 


-6.79 


-7.02 


SS 




-0.44 


-0.45 


-0.44 


SS 


xHqqq 


0.31 


0.32 


0.31 


SS 


x^202 


-1.82 


-1.77 


-1.84 


SD 


a;-D202 


-0.49 


-0.48 


-0.49 


SD 


a; -"202 


0.38 


0.38 


0.38 


SD 


7-1 

a; -'202 


0.023 


0.019 


0.023 


SD 


r)Oii 


1.3 


1.2 


1.3 


SD 


nOll 
x-DlOl 


0.94 


0.90 


0.94 


SD 


trOll 
a; -"101 


-0.66 


-0.62 


-0.66 


SD 


nlll 
a; -^101 


-1.1 


-1.1 


-1.1 


SD 


Rill 

a;-Dl01 


-0.79 


-0.77 


-0.79 


SD 


trill 

a: -"101 


0.55 


0.54 


0.55 


SD 


r)211 


0.29 


0.28 


0.29 


SD 


0211 
a; -"101 


0.19 


0.20 


0.19 


SD 


rr211 
a; -"101 


-0.14 


-0.14 


-0.13 


SD 


a;"-^101 


0.033 


0.036 


0.032 


DD 


rill 

a^-^lOl 


-0.021 


-0.022 


-0.21 


SD,DD 


7-211 
a;-* 101 


-0.023 


-0.024 


-0.025 


DD 
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TABLE V. Trion isovector magnetic moment: (a) AV14, 8 channel (^H | ^H) (b) AV14, 8 
channel ( | ^He) (c) AV14+3BF, 8 channel ( | ^He) (d) AV14+3BF, 22 channel ( | ^He) 



Ref. § 



(a) 



(b) 



(c) 



(d) 



lA 

pair 

pion 

delta 

total 



-2.172 
-0.290 
0.092 
-0.099 
-2.469 



-2.177 
-0.288 
0.096 
-0.100 
-2.468 



-2.174 
-0.285 
0.095 
-0.098 
-2.462 



-2.175 
-0.298 
0.100 
-0.104 
-2.477 



-2.175 
-0.298 
0.101 
-0.105 
-2.477 



TABLE VI. Trion Gamow- Teller matrix element. The entries labeled AHHST are taken from 
Ref. 1^ and use wavefunctions found from the Paris potential. 



Current 


AHHST 


local 


AHHST 


non-local 


AHHST 


AP1-FAP7 


[2a,vert] 


0.000 


0.000 


-0.016 


-0.015 


AP3 


[2a,ret] 


0.000 


0.000 


-0.008 


-0.007 


AP4 


[2a, form] 


0.000 


0.000 


-0.002 


-0.001 


APSPV 


[2a,PS-PV] 


0.000 


0.000 


0.016 


0.015 


AP5 


[2d] 


0.022 


0.019 


0.004 


0.003 


ANPl 


[2c] 


-0.013 


-0.011 


0.000 


-0.001 


ANP2 


[2b] 


0.009 


0.007 


0.003 


0.003 


one-body 


one-body 


0.924 


0.927 


0.000 


0.000 



TABLE VII. Ranges of parameters for the meson exchange currents. 

Parameter range 

3.7 - 6.6 

5piviv/47r 0.70 - 0.95 

Gi 2.2 - 2.6 

A 0-1 

/^jV7v/4vr 0.075 - 0.081 

A^ 1.0 - 1.5 GeV 

/2^^/47r 0.23 - 0.36 



TABLE VIII. Strong form factor fit for two values of rho anomalous coupling ny- The values 
of I^nn/^t^ and g^p^^/Air are 0.081 and 0.95 respectively. The formfactors are given in GeV. 







Ky = 3.7 




Kv = 6.6 


K 


Ap 


Aa 


Ap 


Aa 


1.0 


1.86 


1.09 - 1.86 


1.18 


1.09 - 1.18 


1.2 


2.23 


1.09 - 2.23 


1.25 


1.09 - 1.25 


1.5 


2.76 


1.09 - 2.76 


1.33 


1.09 - 1.33 
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TABLE IX. Results for the rate and spin observables. 





AV i4+ol3r 


\\T-\ /I 1 QT3T7 


AV 14 


TTTJlV/r 
H/r^iVi 




[22] 


[8] 


[8] 


Ref. i 




1502 ± 32 


1498 


1456 


1497 ± 21 


Ay 


0.515 ± 0.005 


0.515 


0.516 


0.524 ± 0.006 


At 


-0.375 ± 0.004 


-0.375 


-0.373 


-0.379 ± 0.004 




-0.110 ± 0.006 


-0.110 


-0.110 


-0.097 ± 0.007 



TABLE X. Results for the rate and spin observables. The one-body currents are taken from 
the 22 channel AV14-I-3BF wavefunctions. 





AV14+3BF 


AV14+3BF 


AV14 




[22] 


[8] 


[8] 


ro [s-i] 


1502 


1501 


1491 


Ay 


0.515 


0.516 


0.518 


At 


-0.375 


-0.375 


-0.374 


A^ 


-0.110 


-0.110 


-0.108 



TABLE XL Contributions of each current (xlO^) to the current amplitudes and effective form 
factors. 



Current 


3a 


,(2) 


Pi 


JQ 


.(1) 

Jx 


Po 


Gy 


Ga 


Gp 


lA 


877 


-326 


10 


15 


200 


820 


835 


1185 


516 


AP1+AP7 


-17 














-17 




AP2+AP6 


-6 


-19 














19 


AP3 


-8 














-8 




AP4 


-3 














-3 




APSPV 


8 














8 




AP5 


22 




1.1 










22 


-1 


ANPl 


-13 




-0.1 










-13 




ANP2 


8 




0.4 










8 




ANP3 


-2 


-5 














5 


ANP4 


77 


4 












75 


-4 


ANP5 


-29 














-29 




ANP6 


2 














2 




VPl 










64 






64 


64 


VP2 










-8 






-8 


-8 


VNPl 










16 






16 


16 


VNP2 










-6 






-6 


-6 


total 


918 


-346 


12 


15 


267 


820 


835 


1300 


601 


EPM 


928 


-372 


11 


15 


241 


839 


854 


1293 


603 



34 



TABLE XII. The split into local and non-local contributions. Only those currents which receive 
contributions from both local and non-local currents are included. 



Current 



>(X103 



ji'^(xl03) 



AP2-FAP6 
AP5 

ANPl 
ANP2 
ANP3 



local 
-4.2 
18.9 

-12.8 
6.2 
-1.2 



non-local 
-2.1 
3.1 
0.2 
2.2 
-0.4 



local 
-12.6 



non-local 
-6.4 



-3.6 



-1.1 



TABLE XIII. Sensitivity of observables O to ^p. 



O 

Ay 

At 
Aa 



gp dO 

O dgp 



„PCAC 

Re 



0.11 
0.37 
0.73 
0.75 



35 



'000 



2.0 



E 
c 

1.0 



0.0 



0.0 



2.0 4.0 
X in fm 



1.0 



'202 



I 0.5 



0.0 L*- 
0.0 



2.0 



X in fm 



4.0 



H 

x"000 



0.0010 



0.0005 



0.0000 



0.0 



2.0 



X in fm 



4.0 



x^lOI 



011 



0.004 



0.002 



0.000 1- 
0.0 



2.0 4.0 

X in fm 




X in fm 



